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Abstract. Combining the second-order entropy flow vector of the causal Israel- 
Stewart theory with the conformal Killing- vector property of Ui/T, where Uj is the four- 
th | velocity of the medium and T its equilibrium temperature, we investigate generalized 

equilibrium states for cosmological fluids with nonconserved particle number. We 
calculate the corresponding equilibrium particle production rate and show that this 
quantity is reduced compared with the results of the previously studied first-order 
| theory. Generalized equilibrium for massive particles turns out to be compatible with a 

""^ ' dependence p oc a~ 2 of the fluid energy density p on the scale factor a of the Robertson- 

Q^. Walker metric and may be regarded as a realization of so-called K-matter. 

PACS numbers: 98.80.Hw, 95.30.Tg, 04.40.Nr, 05.70.Ln 
1. Introduction 

Among the theoretical efforts to understand the physics of the early universe the 
investigation of cosmological nonequilibrium processes is of considerable interest (see, 
e.g., D, |]). A characteristic feature of many nonequilibrium phenomena is the 
nonconservation of particle numbers in interaction processes at high energies within 
the cosmic substratum or due to quantum particle production. Generally, cosmological 
particle creation has been extensively studied in quantum field theory in curved 
spacetime 0, |[] . Progress in nonequilibrium quantum field theory and its cosmological 
applications in the context of fluctuation-dissipation theorems was obtained in || [J 
and references therein. As was realized by Zel'dovich J7], Murphy || and Hu 0, a 
quantum particle production rate may phenomenologically be equivalent to an effective 
viscous pressure. This correspondence has initiated a series of papers discussing a 
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fluid phenomenological approach to cosmological particle production JE)], [II], [T2|, [U| 
[Rj, IT3J, 0, 0, |T|, 0, ID), 0- Special attention has been devoted to "isentropic" 
(or "adiabatic") particle production, characterized by a given, constant equilibrium 
entropy per particle. In such a case the entropy production associated with the creation 
of particles is entirely due to the enlargement of the phase space of the fluid particle 
system and not to fluid dissipative processes which would also enlarge the entropy per 
particle. Particle creation is an intrinsically nonequilibrium process. While a full out- 
of-equilibrium treatment is beyond the present theoretical possibilities, we expect the 



concept of "generalized" equilibrium, previously introduced by one of the authors |22 
to be a step in this direction. "Generalized" equilibrium states share essential features 
with "conventional" equilibrium states. They are, on the other hand, compatible with 
a nonvanishing entropy production. More specifically, this kind of equilibrium relies 
on the existence of a conformal, timelike Killing vector and is compatible with nonzero 
"isentropic" fluid particle production. Despite of the change of the particle number it is 
characterized by a Jiittner distribution function for collisional equilibrium on the level of 
kinetic theory [I 



24, 25 



For massless particles the equilibrium conditions require the 
production rate to vanish and the "generalized" equilibrium coincides with the familiar 
"global" equilibrium of standard relativistic thermodynamics in this limit. Particles with 
nonzero mass for which a "global" equilibrium in the standard sense cannot generally 
exist in an expanding universe, may be in "generalized" equilibrium provided there 
is net production of these particles. The corresponding production rate depends on 
the fluid equations of state and is fixed by the equilibrium properties of the fluid. It 
is highest for nonrelativistic matter in which case it amounts to half the expansion 
rate. In case of spatial homogeneity such kind of equilibrium requires accelerated 
expansion, i.e., it can only be realized in a power-law inflationary universe. This result 
represents a surprising link between equilibrium properties of the cosmological fluid and 
the expansion behaviour of the universe. While the "generalized" equilibrium concept 
is closely related to the properties of a timelike conformal Killing vector it is also based 
on a specific expression for the entropy flow vector. The latter was assumed to be 
of the Eckart type in [22 taking into account first-order deviations from conventional 
equilibrium. It is well known [26, 27, 28, 29, [3(1 that the first-order Eckart theories 



suffer from serious drawbacks concerning causality and stability. If one includes higher- 
order deviations from equilibrium as well, these problems disappear. By now, it is 
generally agreed that any analysis of out-of-equilibrium phenomena in relativity should 
be based on the theories of Miiller f2~6|l , Israel p7fl , Israel and Stewart p8| , or Pavon, Jou 
and Casas- Vazquez []2"9" |, including at least second-order deviations fom (conventional) 
equilibrium. These theories are "causal" in the sense that dissipative signals propagate 
only at subluminal speeds (see |31[] for a recent review). 

It is the aim of this paper to extend the "generalized" equilibrium concept to the 
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second-order theory and to investigate the cosmological dynamics for the case that the 
matter content of the universe is in such kind of equilibrium. We will restrict ourselves 
to (effective) bulk pressures throughout this paper. It should be emphasised that the 
restriction to only one entropy producing phenomenon is more severe in a second-order 
theory than it is in Eckart's theory. In general, the causal evolution equations for the 
viscous pressures and the heat flux are coupled. Consequently, our considerations apply 
to cases where this coupling is not essential. In spatially homogeneous spacetimes, the 
main case of interest here, there is no coupling at all and entropy production is only 
possible due to nonvanishing bulk pressures. It will turn out that the particle production 
rate required by the generalized equilibrium conditions of the second-order theory is 
lower than the corresponding rate in the first-order theory. This strongly reduces the 
possibility of accelerated expansion (power-law inflation) in a spatially homogeneous 
universe as a consequence of the back reaction of equilibrium particle production on 
the cosmological dynamics. The generalized equilibrium conditions of the second-order 
theory may be compatible, however, with a "coasting" cosmology [32], i.e., a universe 
expanding with constant velocity. 

The paper is organized as follows. In section 2 we introduce the concept of 
generalized equilibrium in the context of second-order theories of fluid dynamics. For 
later comparison section 3 briefly recalls the corresponding properties of the first-order 
theory. Section 4.1 discusses the relation between the generalized equilibrium conditions 
and the causal evolution equation for the viscous pressure in the Israel-Stewart theory. 
Section 4.2 demonstrates that the second-order theory admits equilibrium states with 
nonvanishing viscous pressure even for conserved particle number in case Ui/T is a 
conformal Killing vector. In section 4.3 our main results concerning the generalized 
equilibrium production rate for massive particles in the expanding universe are obtained 
and discussed. Section 5 presents the conclusions of the paper. Units have been chosen 
so that c = k B = 1. 



2. The concept of "generalised equilibrium" 

In order to introduce the generalised equilibrium concept, we first recall the standard 
equilibrium properties of relativistic gaseous fluids. A relativistic gas of particles with 
elastic binary interactions, described by Boltzmann's collision integal, is at (local) 
equilibrium if the corresponding one-particle distribution function is of the Juttner 
type. The collision integral vanishes under such circumstances. If, moreover, this 
equilibrium distribution is supposed to obey Boltzmann's equation, the gas is said to 
be at "global" equilibrium. For particles with nonvanishing rest mass this requires the 
existence of a timelike Killing vector, i.e., stationarity. Such kind of equilibrium cannot 
exist in an expanding universe. For massless particles global equilibrium is realised 
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if there exists a conformal timelike Killing vector. This is of cosmological relevance 
since it is compatible with a nonvanishing expansion rate. The concept of "generalised" 
equilibrium is an extension of the global equilibrium concept insofar as it takes into 
account additional interactions within the gaseous fluid which, in general, are not elastic 
and result in a change of the particle number. However, it retains essential features of the 
global equilibrium, namely both the structure of the equilibrium distribution function 
which makes Boltzmann's collision integral vanish (at least within the first-order theory) 
and the conformal Killing vector condition for u^/T. As a consequence, such kind of 
equilibrium becomes possible also for particles with nonvanishing rest mass, provided 
that the particle number increases at a certain rate. In the present paper we shall 
investigate the implications of generalised equilibrium for the second-order theory on 
a phenomenological level. Under these premises we start our considerations with the 
energy-momentum tensor of a fluid with bulk stresses, 

T ik = fnju k + (p + tt) h ik . (1) 

Here, p is the energy density of a fiducial equilibrium state, p is the corresponding 
equilibrium pressure, u % is the fluid four- velocity in the Eckart frame and h lk = g ik + u l u k 
is the spatial projection tensor. The quantity n denotes that part of the scalar pressure 
which is connected with entropy production. The local energy-momentum conservation 
T % \ = implies the balances 

p + Q(p + p + ix) = (2) 

and 

(p + p + tt) ii m + V w (p + tt) = , (3) 

where p = p^u 1 etc. and V m p = /j"p >n . Within the Eckart frame the particle number 
flow vector N % is given by 

N i = nij , (4) 

where n is the particle number density. We will not assume that the fluid particle 
number is conserved. Instead, we admit a source term V in the corresponding balance 
equation which describes the rate of change of the fluid particles: 

N; = nT. (5) 

Inserting here the expression yields 

h + Qn = nT , (6) 

where = u 1 ^ is the fluid expansion. The quantity V describes the net change of the 
particle number either due to internal reactions within the medium |33| or due to particle 



production in strongly varying gravitational fields [0, 0, 0- Usually, a quantity such as 



5 



r is an input quantity in a phenomenological description and has to be calculated from 
the underlying microphysics. Here, it will largely be fixed by the generalized equilibrium 
conditions. 

Limiting ourselves to second-order deviations from equilibrium, the entropy flow 
vector S a takes the form []SSL M \ 

S a = nsu«-^^, (7) 

where s is the entropy per particle and T is the temperature. The quantity r plays 
the role of a relaxation time while ( is identified with the coefficient of bulk viscosity in 
conventional second-order theories. Covariant differentiation of the expression (M) yields 



S* = sN a „ + ns 



From the Gibbs equation (see, e.g., |3 
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one obtains 



Tds = d pd— 

n n 



nTs = —Qtt — (p + p)T 



(8) 



(9) 



(10) 



With the help of Eq. (0) the equation for the entropy production density (||) may be 
written as 



S. a -nTs 



07T p + 1 / 77^ T 
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On a phenomenological level, "generalized equilibrium" is defined as the vanishing of 
S a a — nTs together with the conformal Killing- vector property of Ui/T |22j| : 



si 



nTs 



r > o 



+ 



T 



Ik 



(12) 



The term nTs on the right-hand side of the first relation in Eq. (|T2|) takes into account 
the increase of the entropy of the system due to the enlargement of the phase space. 
Simply due to its very existence each of the created particles contributes an equilibrium 
amount s to the entropy of the fluid system. If there is production of entropy only due 
to the enlargement of the phase space where each particle contributes an equilibrium 
amount s but not due to dissipative processes, the system will be considered to be in 
"generalized" equilibrium. Within the first-order theory the microscopic constituents of 
such a fluid are governed by a distribution function for collisional equilibrium p3| , p^ , [25[ . 
For T = the conditions flT^I ) reduce to the standard equilibrium conditions of relativistic 
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thermodynamics. The conformal Killing- vector (CKV) property of «j/T implies the 
relations 

, i e v a T 



r 3T ' T 

and 

f 



+ Ua = 0, (13) 



(14) 

T 3 V ; 



We point out that the second relation ([13]) implies the vanishing of the heat flux in 
Eckart's theory, equivalent to a vanishing of the corresponding first-order contribution 
to the entropy production within the Israel- Stewart theory. In homogeneous spacetimes 
this relation is trivially fulfilled. 

We will generally assume that the equilibrium variables p and p obey equations of 
state in the form 

p = p(n,T) , p = p(n,T) . (15) 

Differentiating the latter relation and using the balances (|2|) and (|6]) as well as Eq. ([10|) 
we obtain 



Here we have used the general relation 
dp p + p T dp 



(17) 



dn n n dT 

which follows from the fact that the entropy is a state function and the abbreviations 

dp _ (dp/dT) n dp _ ( dp\ 

dp-(dp/dT) n > dT-\dT)j 6C - 

The temperature laws ( p^) and ( Jl6| ) have to be consistent. This requirement fixes 
the production rate. For later comparison we first review the results of the first-order 
theory p2|, equivalent to r = in ([TT|). 



3. The first-order theory 

In the first-order theory the first "generalized" equilibrium condition (^) reduces to (cf 
(0) with r = 0) 

vr = -(p + p)r, (r = 0) , (18) 

equivalent to s = (cf (|10D). For T = this requirement reduces to the condition re = 0, 
dealt with by Bedran and Calvao P5| . 
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Comparing now the temperature laws ( |HD and fll6| ) yields 

r=(l-I^V (r = 0). (19) 



The CKV property of Ui/T fixes the creation rate. An equilibrium according to 
S? a — nsY = is only possible if fluid particles are created at a specific rate given 
by Eq. (|19"D. The production rate T vanishes for p = p/3, i.e. for radiation. This means 
recovering the well-known fact that a fluid obeying an equation of state p = p/3 may be 
in "global equilibrium" if the spacetime admits a timelike CKV. For any other equation 
of state < p < p/3 there exists a specific production rate ( ]T9[ ) which is necessary to 
fulfill the equilibrium conditions. The production of particles at a rate different from 
the rate (|T§| ) disturbs the equilibrium. Except for p = p/3 there is no equilibrium of 
this kind for T = 0. In other words, a fluid with an equation of state in the range 
< p < p/3 may only be in equilibrium if fluid particles are produced at a specific rate 
depending on the equation of state. 

Inserting the relation (|l^) into the condition ([Tj|) fixes ir: 



vr = -(p + p) 1-TTTf , (r=0) . (20) 




The equilibrium conditions naturally imply the existence of an effective bulk pressure, 
i.e., the latter has not to be postulated separately. It is obvious that n as well as T 
vanish for p = p/3. For nonrelativistic matter, however, characterized by equations of 
state p = nT and p = nm + \nT with m ^> T, one obtains 

r = ^e, 7r = -I(p + p)w-| (m»T,r = 0). (21) 

The rate V is half the expansion rate and \tc\ is a remarkable fraction of the energy 
density. With equations of state of the general form (|15D the particle number density n 
and the temperature T are the primary thermodynamic variables of our system. While 
the equilibrium behaviour of the temperature is given by Eq. (|14D independently of the 
equation of state, the particle number density changes according to 

h - = , (r = 0) , (22) 

n 3 op 

where we have combined Eqs. (|6]) and ([05]). Use of ( p~8D and (|T9"| ) in (0) yields 

P = "|(P + P)|, (r = 0). (23) 



For a gas one has [2 
dp 

dp \T J nT V nT 



dp f™) 2 _ 1 + 5 P±P_(e±P)\ (24) 
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Explicit integration of Eqs. (|2"2|) and (|23|) is possible in the limiting cases of radiation 



and nonrelativistic matter. Introducing a length scale a according to 

9 = 3- , (25) 
a 

one recovers the dependences 

n oc a~ 3 , p oc a~ A , (p = p/3 , r = 0) (26) 

for radiation since T as well as tt vanish in this case. For nonrelativistic matter on the 
other hand, one obtains p2[ ^ 



n oc a~ 3/2 , p oc a~ 3/2 , (m > T , r = 0) . (27) 

The temperature behaves according to (|i~4D for any equation of state. Since 
nonrelativistic matter is produced at a rate T = G/2 [cf (|2l|) 1 there are considerable 
deviations from the standard behavior for conserved particle number ( n oc a -3 , T oc a -2 , 
p oc a -3 ). All thermodynamic quantities decrease more slowly than for T = since the 
decay of n, T and p due to the expansion is counteracted by corresponding production 
terms. The differences of the behaviour of n, T and p acoording to the relations ([27]) 
compared with the standard behaviour are consequences of the back reaction of the 
production process on the fluid dynamics. Also the expansion of the universe is modified. 
In a homogeneous and isotropic universe the length scale a coincides with the scale factor 
of the Robert son- Walker metric. Restricting ourselves to the spatially flat case, the scale 
factor obeys the Friedmann equation 

3^ = k P , (28) 

where k is Einstein's gravitational constant. For radiation we recover, of course, a oc t l l 2 . 
But inserting the energy density p of relations ([27]) into we find that the scale factor 
behaves such as 

a oc t 4/3 , (m > T , r = 0) , (29) 

i.e., a > instead of the familiar a oc t 2 ^ 3 with a < for p oc a" 3 , corresponding 
to T = 0. The production of massive particles in equilibrium implies accelerated 
expansion (a > 0), i.e., power-law inflation. In other words, generalized equilibrium 
for nonrelativistic particles is only possible in a power-law inflationary universe. 

4. The second-order theory 

4-1- General relations 



In order to prepare the corresponding investigation of "generalized equilibrium " within 
the second-order theory we first realize [cf (0)] that for T = the first equilibrium 
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condition in (0) is equivalent to 

e 7r + I 7r7 r + ![!l(e + --$--] =0. (30) 

This condition may trivially be fulfilled for tt = 0. However, there exists a different 
possibility, namely a compensation of the terms on the left-hand side of the last equation, 



T7T 



-C9 - -7TT 

s 2 



e + 



(31) 



C T" 
" C " T_ 

We want to point out that this equation for tt does not coincide with the conventional, 
causal evolution equation for the corresponding quantity of the Israel-Stewart theory 
(subscript IS) (see, e.g., 



19, 39 



TT 



(IS) 



1 

— TT T 

2 (is) 



e + 



(32) 



and references therein), which we write down for 
comparison: 

i _ L _ 1 

T ( T 

The equations ([H|) and (|32]) refer to different physical situations. While Eq. (f32|) 
corresponds to an entropy production 

TT 2 

S a , = -f^ , (33) 



which vanishes only for vanishing tt, Eq. (^l|) implies S a a = although it is different 
from zero. The causal evolution equation ([32|) follows from the requirement to guarantee 
S a a > in the simplest possible way with an entropy production given by the expression 
(j3~3|). The viscous pressure 7r. IS) in this case represents deviations from the equilibrium 
S a a = 0. In the present paper we are not interested in such kind of deviations from 
"conventional" equilibrium. Our point is to realize that "generalized" equilibrium within 
the second-order theory does not necessarily imply it = 0. There exists an alternative 
way to make the entropy production vanish, namely a nonvanishing 7r satisfying Eq. 



(pl|). This option relies on a compensation between the contributions to S a a following 
from each of the two terms on the right-hand side of (|7|). One cannot expect, of course, 
that this configuration implies an equilibrium distribution on the level of kinetic theory. 

The bulk pressure n {IS) is a dynamical degree of freedom in Q3"2]), the latter being 
the corresponding evolution equation. Eq. (|31| ) on the other hand, which also may be 
written as 



71 n T\ T 



, (tt + 0) 



(34) 



is an equation for the function r/Q. It will be shown below that the viscous pressure n 
is fixed by the CKV property of «;/T, not yet taken into account in (|3T]), also in the 
second-order theory. 
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Using the generalized equilibrium conditions (|P2|) in ( |TT| ) yields 

97T + {p +p) T + + 2 -^ T -Q = . (35) 



Combination with Eqs. fllQ| ) and ( |25|) allows us to write this condition as 

nTa's = (« 4 y ^) • (36) 

Let us now define a quantity 

1 7r 2 r , , 

s = s , 37 

2nTC, 

with the help of which we may write the entropy flow vector (|7|) as S a = ns*u a . For 
t — > the quantities s* and s coincide. Within the framework of Extended Irreversible 
Thermodynamics (EIT) |37], 38, 3T5| a corresponding expression, 



1 7T T 
S = S — - — , 

(' s > 2nT(' 

is interpreted as nonequilibrium entropy per particle. Since tt, different from tt {IS) , does 
not describe deviations from (generalized) equilibrium (see the discussion following Eq. 
(^)), the quantity ([J7]) is a (generalized) equilibrium variable as well. It may be called 
"second-order generalized equilibrium" entropy per particle. 

With Eq. (g?]) and Ta = const, following from (0) and fl2|), as well as T = N/N, 
where iV = na 3 is the number of particles in a comoving volume a 3 , Eq. (|36| ) may be 
written as 

nTs* = — . (38) 
2 C 

It follows that in the limiting case T = the quantity s* remains constant along the 
fluid flow lines. As we will see below, this does not necessarily imply a vanishing 7r. 

Comparing the rate of change of the temperature ([16]) using s from ([K]) with the 
corresponding rate ( |14D we obtain 

1 dp ( dp\ dp T 



which replaces the expression (|T8j) of the first-order theory. Again, the generalized 
equilibrium conditions establish a relation between the rate V and the bulk pressure 7r. 
While in the first-order theory 7r vanishes with vanishing T, it is obvious from Eq. ([3^) 
that 7r will be nonvanishing in general even for T = 0. Only that part of tt which is 
proportional to T is a "creation" pressure. The first term on the right-hand side of (|39|) 
is independent of the particle production and has to be regarded as a "conventional" 
viscous pressure. This structure of tt is a manifestation of the well-known fact that bulk 
viscosity and particle creation are independent phenomena which, however, may occur 
simultaneously (|T7J . 
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4-2. Particle conservation (T = 0) 



For vanishing T the effective viscous pressure 

dp\ 



reduces to 



7T 




dp 



(40) 



This quantity vanishes for relativistic matter (p = p/3). For any other equation of state, 
however, the equilibrium conditions require a nonvanishing 7r. While Eq. (|35|) may also 
be satisfied with 7r = for T = 0, Eq. (|T4]) may not. Consequently, there exists 
an equilibrium state (vanishing entropy production) within the second-order theory 
which admits a nonvanishing bulk pressure 7r. This apparently new result relies on 
the compensation effects mentioned below 



50p and in between (^) and (|34|). We recall 
that this case will certainly not be realised microscopically by a distribution function for 
collisional equilibrium. The part of 7r which is not due to T represents a "conventional" 
bulk viscous pressure. These considerations show that within the second-order theory 
the set of equilibrium relations ( |T2D admits qualitatively new configurations compared 
with the first-order theory. In fact, the case T = discussed here is a consequence of the 
standard equilibrium requirement = under the condition that u a /T be a conformal 
Killing vector. 

For nonrelativistic matter (p = nT, p = nm + |nT, m ^> T) one obtains 



7T 



— nT 
2 



i.e. 



|tt| = p/2, although Sf a = 0. 
[cf (|38|) 1 is equivalent to 



12^ 

a 



Using the expression (|4i~| ) for n in 




(41) 

the property 
(42) 



In 



(43) 



Integration yields 

r 12 

( nT \ao, 

The quantity t/( is only determined up to a constant. We have used here this freedom 
to make t/( vanish for some initial value a = ao- Since t/( vanishes identically in 
the first-order theory this may be interpreted as "switching on" the second order at 
a = oq. The function r/C then describes deviations from the corresponding first-order 
description. In the Israel-Stewart theory the function r/C is known to be related to the 
viscous part q, of the sound velocity in a viscous fluid by 

1 



40 



C 



(44) 
p + pr 

This expression follows also in the appropriate limit from formula (59) of [f|T| with the 
identifications nF — > p+p and (3q = It is remarkable that the system of perturbation 
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equations (53) - (56) in from which the velocity ( [P]) is obtained, does not change if 
Eq. fl32"D of the general analysis is replaced by Eq. (|31|). Therefore, the expression ([14]) 
may be used in the present context as well. We associate q, with the effective viscous 
part of a propagation velocity in a fluid in generalized equilibrium. The quantity c\ 
diverges for r — > (vanishing second order) reflecting the "acausal" nature of the first- 
order theories. The more familiar adiabatic part of the sound velocity c s is given by 



2 / w dp T {dp/dTf 

s {a ~' ad p + pdn p + p dp/dT ' 1 ' 




The sound in a viscous medium propagates with a subluminal velocity v where |40[ 

v 2 = c 2 s + c 2 b <l. (46) 
With the equations of state for nonrelativistic matter we find 

c; 2 = |^, c- = i2^mm. (47) 



5T ' T Va , 

For a sufficiently larger than ao (recall that a = ao represents the (unphysical) first-order 
limiting case) we have ^ < 1, t| < 1 and < c 2 s1 i.e., the "adiabatic" part exceeds 
the viscous one which is consistent with \n\ < p according to (fy]). This completes 
our discussion of equilibrium states of the second-order theory for the case of vanishing 
particle production. 



4-3. Particle production (T > 0) 



According to fl38|) the generalized entropy per particle s* does not remain constant for 
T > 0. Taking into account Eq. ([37]) and iV = na 3 , Eq. ( ^8|) may be written as 



[N {**-*))' 



-Ns 



(Ah 



The first-order case corresponds to s* = s, implying s = 0. We will integrate Eq. (fi"8|) 
for the limiting cases of pure radiation and nonrelativistic matter. In the first case 
(p = nT, p = 3nT) the expression (|39| ) for the viscous pressure reduces to 



7T 



~ti~~ = —3nT- 



(49) 



dnQ 

Use of the last equation together with relation fllTD and p = nT in (|T0|) yields 

Ns = -NT = -N (50) 
for the right-hand side of ( |8 ) which, consequently, may be written as 

1 7T 2 r N 



AM 1 + 



2nTC, 







(51) 
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with 7r from (f|9). Assuming again the funtion r/C to vanish for some initial time, i.e. 
( T /C)o = 0; the solution of (|5T|) is 



In 2 r N-N, 



(52) 



2nT( N 

where N is the initial particle number. For N > N the last relation may only be 
fulfilled if the ratio t/( is negative. Since formula (|44|) necessarily implies a nonnegative 
value of this ratio, we conclude that the only physical solution of (|5lD is it = T = 0. 
This confirms the corresponding result of the first-order theory according to which the 
production of relativistic particles is incompatible with the conditions of generalized 
equilibrium. 

For nonrelativistic matter (jp = nT, p = nm + |nT m ^> T) the viscous pressure 
fl3"9|) becomes 

vr = -~nT - p^ , (m » T) . (53) 

For T = we recover the limiting case (|4TD. Introducing the expression (|53| ) into ([10|) 
yields 

s = -6 - T , (m > T) . (54) 

The condition s = correctly reproduces the relation T = 0/2 [cf fl2l"D l of the first-order 
theory. For Ns on the right-hand side of (^) we obtain 

Q 

Ns = -N— — N , (m > T) . (55) 

2 a 

To integrate this expression we assume T to have the form 

T = ^9 , (m » T) , (56) 

with a constant parameter f3. The case (3 = 1 corresponds to the first-order theory. 
Differences between the second- and first-order theories are then characterized by 
deviations from (3 = 1. With the ansatz (|56| ) integration of the particle number balance 
(|) yields 

3/3/2 



N = No 



— ) (m>T) . (57) 

do J 



Using this law for iV in Eq. (^BJ), integration of the latter results in 



* Nsdt = ^—^ (N - N ) , (m > T) . (58) 

t (3 

Applying now the definition (|3?D of s* and assuming again r/C to vanish initially (see 
the discussion following (|4"3"D), the integral of (^) turns out to be 

jr^r _ 1-PN-Nq 

2nT( (3 N • 1 J 
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The limit (3=1, i.e. F = 0/2, corresponds to r = of the first-order theory. The 
case (3 > 1 is unphysical since it is incompatible with a positive t/(. Consequently, 
the second-order theory requires (3 < 1. In other words, the generalized equilibrium 
conditions of the second-order theory imply a production rate which is less than the 
corresponding rate of the first-order theory. To be more specific, we assume (3 to be 
smaller but still of the order of 1. In such a case ir is approximately given by 

7r — ^nm , (m^$> T) , (60) 

i.e., it is entirely determined by the particle production. Use of r/C = [(p + p) c 2 ] 1 ~ 
[rime 2 ] 1 [cf (f44])1 and relation ( |60"D allows us to write the left-hand side of (|59| ) as 

1 7T 2 T B 2 m , , , 

2^? = k T ' {m >y T) ' (61) 

Again, the case r — > corresponds to the (unphysical) behaviour c 2 — > oo. Combining 
Eqs. ( |5D| ) and ( |5I| ) and assuming iV 3> A^ , i.e. considerable particle production, we 
obtain 

, 1 3 3 m m . . , 



Recalling from Eq. (|47|) that the adiabatic part of the sound velocity is characterized 
by c 2 1 for T m, the condition f 2 < 1, restricting the sound propagation to 
subluminal speeds, reduces to 



m 8(1-/3) 

c '<! =► t < ^^- (63) 



The latter inequality has to be fulfilled for m/T ^> 1. These simultaneous requirements 
restrict the possible /3-values and, according to (|5E| ) limit the production rate. While 
the generalized equilibrium conditions of the first-order theory imply the result (3 = 1 
for nonrelativistic particles, (3 is not uniquely fixed within the second-order theory. The 
conditions m ^> T and QBlf ) determine only a range of possible /^-values. With (3 = 1/2, 
e.g., we find 8 (1 — (3) //? 3 = 32, a value perfectly admissible for m/T. For /3 = 2/3 one 
has m/T < 9 which means that m still may be almost an order of magnitude larger 
than T. For (3 = 3/4 we obtain m/T < 4.74, a value already hardly compatible with 
m ^> T. For values of (3 still closer to 1 the upper limit for m/T becomes even lower. 
We conclude that the causality requirement of the second-order theory restricts (3 to 
values of the order of 2/3, i.e. 

2 6 1 

(3<- T<— \n\<-nm, (m > T) . (64) 

3 3 3 

The limiting value T = 0/3 corresponds to n = —nm/3. Use of the latter in the energy 
balance @ implies 

p oc a~ 2 , (m>T) . (65) 
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The second-order theory predicts p to decay faster than the first-order theory [cf (pTf)]. 
A behaviour of the energy density such as given by Eq. (|65|) is characteristic for so-called 
"K-matter" . Consequently, generalized equilibrium for massive particles within the 
second-order theory may be regarded as realization of K-matter. In a spatially flat 
universe the behaviour ( |Bo| ) implies a oc t, i.e. a = (coasting cosmology). While for 
massive particles the generalized equilibrium of the first-order theory is connected with 
power-law inflation a oc a 4//3 , i.e. accelerated expansion, the second-order theory reduces 
the acceleration of the expansion. The first-order theory uniquely fixes the power q in 
a oc t q for the scale factor a of a spatially flat universe of massive particles to q = 4/3. 
The second-order theory restricts this power to q < 1. In other words, the first-order 
theory overestimates the generalized equilibrium value for the particle production and 
the back reaction of the latter on the cosmological dynamics. 

5. Summary 

We have defined generalized equilibrium in second-order thermodynamics by the set 
of conditions (|i~2"|) together with the expression ([F[). For p = p/3 (radiation) this kind 
of equilibrium reduces to the standard "global " equilibrium of relativistic thermody- 
namics and requires particle number conservation both in the first- and second-order 
theories. For any other equation of state the first- and second-order generalized equilib- 
rium states are different. While the first-order equilibrium necessarily implies nonzero 
particle production (except for p = p/3) the second-order theory admits (generalized) 
equilibrium both for vanishing and nonvanishing production rates. In the first case (no 
particle production) vanishing entropy production is compatible with a nonvanishing 
viscous pressure. In the second case (nonzero particle production) the equilibrium par- 
ticle production rate is lower than the corresponding first-order quantity. Consequently, 
the back reaction of the production process on the cosmological dynamics is weaker than 
predicted on the basis of the first-order theory. Especially, the property of power law 
inflation a oc in a spatially homogeneous and flat universe filled with nonrelativistic 
matter in first-order generalized equilibrium does no longer hold in the second-order 
theory. The maximally possible production rate in the second-order theory is of the or- 
der of the Hubble expansion rate if = 0/3 compared with 0/2 of the first-order theory. 
This limiting value (r « H for nonrelativistic matter) of the second-order theory implies 
the dependence p oc a~ 2 (a = const) for the energy density, a behavior characteristic for 
K-matter. A universe of massive particles in generalized second-order equilibrium may 
realize a "coasting" cosmology. 
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